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Abstract 

The Hawking emissivities for the scalar-, vector-, and tensor- mode bulk gravi- 
tons are computed in the full range of the graviton's energy by adopting the 
analytic continuation numerically when the spacetime background is (4 + Tri- 
dimensional non-rotating black hole. The total emissivity for the gravitons is 
only 5.16% of that for the spin-0 field when there is no extra dimension. How- 
ever, this ratio factor increases rapidly when the extra dimensions exist. For 
example, this factor becomes 147.7%, 595.2% and 3496% when the number 
of extra dimensions is 1, 2 and 6, respectively. This fact indicates that the 
Hawking radiation for the graviton modes becomes more and more significant 
and dominant with increasing the number of extra dimensions. 



* Email : dkpar k@hep . kyungnam . ac . kr 



1 



Recent quantum gravity such as string theories [1] or brane-world scenario [2] generally 
requires the extra dimensions to reconcile general relativity with quantum physics. Espe- 
cially the modern brane-world scenarios predict the emergence of the TeV-scale gravity, 
which opens the possibility to make tiny black holes by high-energy scattering in the fu- 
ture colliders [3]. In this reason much attention is paid recently to the effect of the extra 
dimensions in the black hole physics. 

The absorption problem and Hawking radiation for the spin 0, 1/2 and 1 particles in 
the background of the (4 + n)-dimensional Schwarzschild black hole were explored in Ref. 
[4]. The numerical calculation supports the fact that the black holes radiate mainly on the 
brane. In fact this was pointed out by Emparan, Horowitz and Myers(EHM) in Ref. [5] by 
making use of the higher-dimensional black body radiation. This claim was also supported 
in the background of the higher-dimensional charged black hole [6]. 

More recently, this issue was re-examined when the situation is different. If, for example, 
black hole has a rotation, there is an important factor we should consider carefully called 
superradiance [7] . The superradiance in the background of the higher-dimensional black hole 
was examined for the bulk fields [8] and brane fields [9]. However, numerical calculation has 
shown that in spite of the consideration of the superradiance EHM claim still holds due to 
the incrediably large difference in the energy amplification for the brane field and bulk field 
[10]. 

There is an another factor we have not considered thoroughly. This is an Hawking radi- 
ation for the higher-spin particles like graviton. Since the graviton is not generally localized 
on the brane unlike the usual standard model particles, the argument of EHM should be 
carefully re-checked in the graviton emission. Generalizing the Regge- Wheeler method [11], 
the various gravitational perturbations were studied in Ref. [12] in the background of the 
higher-dimensional Schwarzschild black hole. Using the radial equations derived in Ref. 
[12], the low-energy and high-energy behaviors for the bulk graviton absorption and emis- 
sion spectra were recently studied [13]. The graviton emission on the brane is also examined 
using an axial perturbation [14]. In Ref. [14] it was argued that the graviton emission can 
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be dominant one in the Hawking radiation when there are many extra dimensions. We 
would like to explore this issue again in the bulk emission. In the following we will compute 
the absorption and emission spectra for the various graviton modes numerically. We will 
show the emission rates are generally enhanced when the number of extra dimensions, say 
n, increases. However, the increasing rate for the gravitons is much larger than that for the 
spin-0 field. For example, the total emissivities for the gravitons is only 5.16% of that for the 
spin-0 field in four dimensions. However, this ratio becomes 595.2% when n = 2 and 3496% 
when n — 6. Thus the Hawking radiation for the higher-spin fields becomes dominant when 
the extra dimensions exist. 

We start with (4 + n)-dimensional Schwarzschild spacetime whose metric is 



ds 2 = -hdt 2 + h~ l dr 2 + r 2 dil 2 



n+2 



(1) 



where h — 1 — (rf//r) n+1 and the angle part dQ^^ is a spherically symmetric line element 
in a form 



dtt 2 n+2 = d0 2 n+l + sin 2 9 n+1 



d6 2 n + sin2 °n ( ■ ■ ■ + sin 2 6 2 (d# 2 + sin 2 My? 2 ) 



(2) 



It is well-known [12] that when n^O, there are three gravitational metric perturbations, 
i.e. scalar, vector and tensor perturbations. For the vector and tensor modes the radial 
equation reduces to the following Schrodinger-like expression 



(3) 



where r* is a "tortoise" coordinate defined dr/dr* = h and the effective potential Vvt is 



h 



VT 



+ n+ 1) + 



n{n + 2) k{n + 2) 2 (r H \ n+l 



(4) 



k 



4 4 - r 

— 1 for tensor mode 
3 for vector mode 
with £ > 2. 

The radial equation for the scalar mode also reduces to the Schrodinger-like expression. 
However, the effective potential Vg is comparatively complicate as following 



h qy 3 + py 2 + wy + z 
rM[2m + (n + 2){n + 3)yf 



(5) 



where y = 1 — h = (r H /r) n+1 and 



(6) 



m = £(£ + n + l) - (n + 2) q = (n + 2) 4 (n + 3) 2 

z = 16m 3 + 4m 2 (n + 2)(n + 4) 

p = {n + 2)(n + 3) [4m(2n 2 + 5n + 6) + n(n + 2)(n + 3)(n - 2) 
w = -12m(n + 2) [m(n - 2) + n(n + 2)(n + 3)] 

with £ > 2. 

In the 4d limit the vector mode corresponds to the gravitational axial perturbation [15], 
whose effective potential reduces to 



Mr) = y 2 



(7) 



and the scalar mode corresponds to the gravitational polar perturbation [16] with 



Vp(r) 



h 



8A 2 (A + 1) + 12A 2 (^)+18A (^-J ^ + 9 



(8) 



(2Ar + 3r H ) 2 

where A = (£ — 1)(£+2)/2. There is no correspondence of the tensor mode in four dimension. 
The most striking result in the 4d gravitational metric perturbations is the fact that the 
effective potentials Va and V P) which look completely different from each other, are expressed 
as a single equation [17] 



v P , A (r) = ±£ + ff 2 + K f 



(9) 



where (3 = 3r H , k = 4A(A + 1) and / = h/r(2\r + 3r H ). In fact, this relation was found when 
Newman-Penrose formalism [18] is applied to the 4d gravitational perturbations. Making 
use of this explicit relation, one can show that the effective potentials Va and Vp have the 
same transmission coefficient. This fact also indicates that the absorption and emission 
spectra for the vector mode graviton and scalar mode graviton by a Schwarzschild black 
hole are exactly same in four dimensions. 



Now, we would like to discuss how the absorption and emission spectra are computed. We 
first consider the vector and tensor modes of the bulk graviton. Defining the dimensionless 
parameters x = our and xr = urn, one can rewrite Eq.(3) in the following 

4a; 2 (x n+1 - ^ +1 ) 2 TT + 4 ( n + 1 ) x h + ' x ~ x h +1 ) ^ (10) 

+ [4x 2n+4 - (x n+l - x n H +1 ) [(21 + n){2£ + n + 2)x n+l - k{n + 2) 2 x n + 1 )] * = 0. 

One can easily show that if ^ is a solution of Eq.(10), its complex conjugate is also 
solution. Eq.(10) also guarantees the Wronskian between them is 

dfy d®* Cx n+1 
[ ' lx ~ dx dx x^ - x n H +1 1 } 

where C is an integration constant. 

Now, we would like to consider the solution of Eq.(10), which is convergent in the near- 
horizon regime, i.e. x ~ xh- Since x = xh is a regular singular point, this solution can be 
derived by a series expression as following 

oo 

M^%) = e Ar,lnMfll E^-^) Ar - (12) 



7V=0 

Inserting Eq.(12) into (10) easily yields 

A„ = -^. (13) 

The sign of A n is chosen for the near-horizon solution to be ingoing into the black hole. 
The recursion relation for the coefficients de,N can be obtained too. Since it is ,of course, 
n-dependent and lengthy, we will not present it here. The Wronskian (11) enables us to 
derive 

where g n/ = d efi - 

Next, we would like to discuss the solution of Eq.(10), which is convergent in the asymp- 
totic regime. This is also expressed as a series form by inverse power of x as following: 



r nA±) {x, x H ) = {±i)^xe^{x - x H )^ ]T r N (±)X^ N+1) (15) 

N=0 

with r (±) = 1. The solutions J- n ,t{+) and J- n ,t(-) are respectively the ingoing and outgoing 
solutions. Using Eq.(ll), one can show easily 

2ix n+1 

W[F nA+) , F^-^ = — — — -— T . (16) 

X X jj 

The recursion relation for the coefficients T/v(±) is straightforwardly derived, which is too 
lengthy to present here. 

Now, we would like to discuss how the coefficient g n ,t{= d>e,o) is related to the scattering 
amplitude. Since the real scattering solution, say 7Z n / should be ingoing wave in the near- 
horizon regime and the mixture of ingoing and outgoing waves at the asymptotic regime, 
one may express it in the form: 

n n/ X ~^ H g n/ (x - x H ) x " [1 + 0(x - x H )\ (17) 

n n/ 7^ r v ^-J Q ^ 

-ix+\„\n\x-x H \ _ ^_iy+% g £( y X}j)e ix ~ Xn ^ ll \ x ~ XH \ + O (— 



X 

\x, 

where Sn^xn) is a scattering amplitude and Q n / is a quantity related to the multiplicities 
M n ,i as following: 
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Q n e = _ -J\f £ . (18) 

W (n + 4)(n + l) ' V ; 

The multiplicities A/^?, and Afffl for the scalar-, vector- and tensor-mode bulk gravi- 

tons are given by [13] 

(5) _ (2£ + n+l)(£ + n)\ 

N ^ ~ (19) 

f{y) _ £(£ + n+l)(2£ + n+l)(£ + n-l)\ 

n ' e ~ (£+l)\n\ 

f{T) _ n(n + 3)(£ + n + 2)(£-l)(2£ + n+l)(£ + n- 1)! 

n ' e ~ 2(£ + l)!(n+l)! 

Introducing a phase shift <5 n ^ as S n j = e 2tSn - e , one can rewrite the second equation of 
Eq.(17) in the form 



K n , t ^=T (1±?) Q n/ e^ (20) 



x sin 



x + i\ n In \x - x H \ - ^ \£ + ^ ) <)„ 



+ 



The near-horizon behavior of the real scattering solution, i.e. the first one of Eq.(17), 
implies that the Wronskian W\TZ* n plZ n ^ x is same with Eq.(14). However, Eq.(20) implies 

9™ / 1 _|_ »j \ T n+1 

[K p n n/ ] x = -*-r 2 (-I-) ^ e -^smh2^— (21) 

where S n> £ = rj n ^ + if3 n ^. Equalizing these two Wronskians yields a relation between g n ^ and 
P n / as following: 

\9nA 2 = (i±^) Q 2 , (l - e" 4 ^) . (22) 

Thus one can calculate the transimission coefficient 1 — l-Sy^ 2 = 1 — e _4/3n -* if g n j is known. 

Now, we would like to explain how to compute g n ^. For the explanation it is convenient 
to introduce a new radial solution lZ n ^(x), which differs from TZ n ^{x) in its normalization 
in such a way that 

n n/ (x, x H ) * ~* h (x - x H ) x " [1 + 0(x - x H )\ . (23) 

Since ^F n ,e{±) are t wo linearly independent solutions of Eq.(10), one can express 1Z n j as a 
linear combination of them as following 

£n,/(a;, xjy) = f$ (x H )Fn4{+) (x, x H ) + fiy (xH)Fn,e(-) ( x > x h) (24) 

where are called jost functions. Using Eq.(16) one can compute the jost functions as 
following 

= ± 2 . xn+1 H w[f nA±) ,n n , e ] x . (25) 

Inserting the explicit expressions of T n ,t{±) in Eq.(15) into Eq.(24) and comparing it with 
the second equation of Eq.(17), one can derive the following two relations 
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SuA x h) = -py- — - (26) 
fnA x n) 

o — — 1 1 

Combining Eq.(22) and (26), one can express the greybody factor (or transmission coefficient) 
of the black hole in terms of the jost function in the following 

1 - i^i 2 = rFv (27) 

\Jn,e I 

The partial absorption cross sections for the vector and tensor graviton modes by the 
(4 + n)-dimensional Schwarzschild black hole are given by 



V 2 J Qn > e r n + 2 lf (-) l2 - ( 2§ ) 
X H \Jn,e I 

Applying the Hawking formula [19], one can compute the bulk emission rate, i.e. the energy 
emitted to the bulk per unit time and unit energy interval, as following 



r 

dui 



2 „ + 3 x( „ + 3> /2r + (n + 4)(n + 



(29) 



where a a b s is a total absorption cross section and Th is an Hawking temperature given by 
Th — (n+ l)/47rrH- Thus one can compute the absorption and emission spectra in the full 
range of uj if the jost functions are computed. 

One can compute the jost functions by matching the near-horizon and asymptotic solu- 
tions. This is achieved by the analytic continuation [6,10] using a solution of Eq.(10) which 
is convergent at the arbitrary point x = b. This solution is also expressed as a series form 
as following 

CO 

<p n ,e( x > x h) = ( x ~ x n) Xn E D n( x - b) N . (30) 

N=0 

The recursion for the coefficients Dn can be obtained by inserting Eq.(30) into (10). Since 
it is too lengthy, we will not present it here. Using a solution (30) one can increase the 
convergent region for the asymptotic solution from the near-horizon regime and decrease for 
the asymptotic solution from the asymptotic regime. Repeating the procedure eventually 
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makes the two solutions which have common convergent region. Then one can compute the 
jost functions by making use of these two solutions and Eq.(25). 

Now, we would like to discuss the case of the scalar mode graviton. Introducing the 
dimensionless parameters x = uir and xh = ojth again, one can rewrite the radial equation 
for the scalar mode as 

Ax 2 (x n+1 - x n H +1 f [2mx n+1 + (n + 2) (n + 3)x n H +1 ] 2 ^ (31) 
+4(n + l)x n H +1 x (x n+1 - x n H +1 ) [2mx n+1 + (n + 2)(n + 3)x" +l1 '"' " * 



dx 



Ax 



2n+4 



2mx n+1 + (n + 2)(n + 3)x n H +1 



- (x n+1 - x r ^ 1 ) {zx 3n+3 + wx n H +1 x 2n+2 + px% +2 x n+1 + gx| i+3 } 



* = 0. 



The calculational procedure is exactly same with the case of vector or tensor mode. All 
Wronskians in Eq.(ll), (14) and (16) as well as the multiplication factor in Eq.(13) are 
exactly same with the previous case. The differences arise only in the recursion relations for 
the coefficients de t N, t n(±) and and the multiplicity defined in Eq.(19). Although the 
recursion relations in the scalar mode case are much more complicate and lengthy, they do 
not make any difficulty in the numerical calculation. In actual numerical computation we 
summed up the partial modes for 2 < I < 12 to plot the absorption and emission spectra of 
the graviton modes. 

Fig.l shows the total absorption and emission rates for the vector (or axial)-mode graviton 
and scalar(or polar)-mode graviton when n — 0. To compare the graviton spectra with those 
for the low-spin field, we plot the absorption and emission rates for the spin-0 scalar field 
together. As Eq.(9) implies, the absorption and emission spectra for both graviton modes 
are exactly identical. As Fig. 1(b) indicates, the emission rates for the gravitons are much 
smaller than that for the spin-0 field. The total emission rates for both gravitons are only 
5.16% of that for the spin-0 field. Thus, in Ad the Hawking radiation for the low-spin fields 
such as standard model fields is dominant one. However, the situation is drastically changed 
when the extra dimensions exist. 
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Fig. 1(a) 
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scalar field 

scalar-mode graviton 
& 

\ vector-mode graviton 



2 
Fig.1(b) 

FIG. 1. The absorption (Fig. 1(a)) and emission (Fig. 1(b)) spectra for the scalar- and vec- 
tor-mode gravitons when n = 0. For a comparision the spectra for the spin-0 field are plotted 
together. As Eq.(9) implies, the spectra for both gravitons are exactly identical. Fig. 1(b) indicates 
that the emission rates for the gravitons are negligible compared to that for the spin-0 field when 
there is no extra dimension. 



Fig. 2 shows that the total emission rates for the scalar-, vector-, and tensor-mode gravi- 
tons when n = 1, 2, and 6. Like Fig.l the emission spectrum of the spin-0 field is plotted 
together for a comparision. As Fig. 2 indicates, the emissivities for all fields generally in- 
crease with increasing n. However, the increasing rates for the gravitons in the Hawking 
emissivities with increasing n is much larger than that for the spin-0 field. 
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spin-0 


scalar-mode / spin-0 


vector-mode / spin-0 


tensor-mode / spin-0 


Ad 


0.000297531 


0.0258 


0.0258 




5d 


0.00104968 


0.8115 


0.60405 


0.06161 


6d 


0.002679 


3.02476 


2.37418 


0.55351 


lOd 


0.118857 


4.94863 


13.4278 


16.5793 



Table I 



Table 1 shows the relative emissivities when n — 0, 1, 2 and 6. When n — 2, the emission 
rate for the scalar mode graviton is almost three times than that for the spin-0 field. When 
n — 6, the emission rate for the tensor-mode graviton is seventeen times than that for the 
spin-0 field!!! The remarkable fact is that the ratio of the tensor-mode graviton to the spin-0 
field increases rapidly with increasing n compared to other graviton modes. 

In this letter we computed the absorption and emission spectra for the various modes of 
the bulk gravitons in the higher-dimensional non-rotating black hole background. The total 
emissivities for the gravitons are only 5.16% compared to that for the spin-0 field in four 
dimensions. However, this ratio increases rapidly when the extra dimensions exist. When, 
for example, n — 1, 2 and 6, this ratio goes to 147.7%, 595.2% and 3496%, respectively. 
This fact indicates that the emission of the higher-spin field like gravitons becomes more and 
more dominant and has a experimental significance when there are many extra dimensions. 
Same conclusion was derived in the emission on the brane [14] . Thus, it is important in our 
opinion to re-examine the EHM argument 11 black holes radiate mainly on the brane" with 
considering the Hawking radiation for the gravitons. We hope to study this issue in the 
future. 
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FIG. 2. The emission spectra for the scalar-, vector- and tensor-mode gravitons when n = 1 
(Fig. 2(a)), n = 2 (Fig. 2(b)) and n = 6 (Fig. 2(c)). The spectra for the spin-0 field are plotted 
together for a comparision. In general the emissivities for all fields increase with increasing n. 
However, the increasing rate for the gravitons are much larger than that for the spin-0 field. Thus, 
Hawking radiation for the gravitons becomes more and more significant when the extra dimensions 
exist. 
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